Critical exponent in the magnetization curve of quantum spin chains 
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The ground state magnetization curve around the critical magnetic field H c of quantum spin chains 
with the spin gap is investigated. We propose a size scaling method to estimate the critical exponent 
5 defined as m ~ \H — H^ 1 ^ 5 from finite cluster calculation. The applications of the method to the 
S — 1 antiferromagnetic chain and S = 1/2 bond alternating chain lead to a common conclusion 
8 = 2. The same result is derived for both edges of the magnetization plateau of the 5* = 3/2 
antiferromagnetic chain with the single ion anisotropy. 
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The magnetization curve of quantum spin chains shows 
various nontrivial behaviors due to quantum effects. In 
the spin-gap systems, where a finite energy gap exists 
in the spin excitation spectrum, the gap is controlled 
by applied magnetic field H through the Zeeman term 
in the Hamiltonian. The typical examples are the an- 
tiferromagnetic chain with integer spin called Haldane 
magnets jlj], spin Peierls systems and spin ladders etc. 
In these systems a phase transition occurs at the criti- 
cal field H c corresponding to the amplitude of the gap 
the system has the nonmagnetic ground state and 
a finite gap for H < H c , while the magnetic ground 
state and no gap for H > H c . The transition was 
observed in the magnetization measurements on some 
quasi-one-dimensional materials; for example, an S = 1 
antiferromagnet Ni(C 2 H 8 N 2 ) 2 N02(C104) |,|, abbrevi- 
ated NENP, and a spin-Peierls compound CuGeC>3. [7]]|] 

In our previous work Q we presented a method to 
derive the ground-state magnetization curve in the ther- 
modynamic limit from the finite-cluster calculation by 
the size scaling based on the conformal invariance. || 
The obtained curve of the 5 — 1 antiferromagnetic chain 
successfully realized the experimental results of the mag- 
netization measurements on NENP qualitatively. The 
method was also applied to get theoretical magnetiza- 
tion curves of some other one-dimensional spin systems. 
p0| , p| However, the critical behavior near H c cannot be 
investigated by this method, because it can yield too few 
points near H c to determine the critical exponent of the 
magnetization curve by the standard curve fitting. 

In general, except for the Kosterlitz-Thouless transi- 
tion p^ |, the magnetization m near the critical field be- 
haves like 



(1) 



for the second-order phase transition. The critical ex- 
ponent 6 is an important quantity to determine the uni- 
versality class of the phase transition which does not de- 
pend on any detailed properties of each system. For the 
5 = 1 antiferromagnetic chain the exponent was deduced 
as 6 = 2 from some effective Hamiltonian theories. [^3|Jl4| 



For the 5 = 1/2 bond alternating chain the bosonization 
method gave the same result. jl5| 6 = 2 was also de- 
rived from the fermionic excitation with the dispersion 
k 2 which was numerically verified to be a good picture 
for both systems. |jj.^J In addition the argument of the 
equivalence between the magnetization process of anti- 
ferromagnetic chains and some integrable models of the 
crystal-shape profile lead to the same conclusion. |l6| In 
any theories giving 6 = 2, however, the original spin 
Hamiltonians were mapped into other solvable models 
with some crucial approximations. Thus it would be im- 
portant to estimate 6 for the original Hamiltonian di- 
rectly in some numerical ways, to test these effective the- 
ories and to investigate unknown systems. 

In this paper we propose a size scaling method to es- 
timate the critical exponent 6 of quantum spin chains 
using the result of the finite cluster calculation. In or- 
der to examine the validity of the method, we apply it 
to the 5=1 antiferromagnetic chain and the 5=1/2 
bond alternating chains. In addition a recent topic on the 
magnetization plateau of the 5 = 3/2 antiferromagnetic 
chain with anisotropy is investigated by the method. 

At first we consider the 5=1 antiferromagnetic 
Heisenberg chain for the explanation of the method. The 
following argument is easily to be applied to more gener- 
alized models. To investigate the magnetization process 
we consider the Hamiltonian 



H = Ho + Hz, 

L 

H =*^2 S i ' S i+i> 

i 

L 

Hz= — Sj, 



(2) 



under the periodic boundary condition. We restrict us on 
even-site systems to avoid the frustration. Throughout 
we use the unit such that g/iB = 1- For L-site systems, 
the lowest energy of Ho in the subspace where J2j Sj = 
M (the macroscopic magnetization is m = M/L) is de- 
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noted as E(L,M). We assume the asymptotic form of 
the size dependence of the energy as 



±.E(L,M)~ e (m)+C(m)^ s (L 



(3) 



where e(m) is the bulk energy and the second term de- 
scribes the leading size correction. We also assume that 
C(m) is an analytic function of m. For gapless cases 
the conformal field theory predicted 9 — 2. || Since the 
method works better for faster convergence of the size 
correction as shown later, we can also accept the expo- 
nential decay like e~ L ^ which is reasonably expected for 
the ground state of the spin gap systems instead of l/L 6 . 
We neglect the m-dependence of 9 because it gives only 
higher order corrections which does not change the main 
result. If the bulk system has the critical behavior de- 
scribed by the form ([!]), m-dependence of the energy e(m) 
near to = should have the form 



e(m) - e(0) + H c m + Am 



s+i 



(4) 



where A is a positive constant and we assume S > 0. Now 
we put M =0, 1 and 2 into the form (|J) and use (||). If L 
is sufficiently large, C(m) can be expanded with respect 
to m as C(l/L) ~ C(0) + C'(0)1/L + 1/2C"(0)1/L 2 • • • 
for M = 1. Thus we get the forms 



E(L,0) 
E(L,l) 

E(L,2) 



Le(0) + C(0)- 



1 



£8-1 > 

Le(0) + H C + A-^ 



Le(0) + 2H C + 2 5+1 Ajj 
+C(0)-^ T + 2C'(0)- 1 



2C"(0) 



L e+i 



If we define the quantity 

f(L) = E(L,2)+E(L,0)-2E(L,l), 
the asymptotic behavior of f(L) becomes 

f(L)~A(2 s+1 - 2)-L + C"(0)-^ T (L-oo). (6) 



(5) 



When the second term of (|J) converges faster than the 
first one, the exponent 5 can be estimated from the size 
dependence of f(L). Thus the necessary condition under 
which the method gives the correct value of 5 is 



><5-l. 



(7) 



Therefore we have to check the condition that E(L, 0)/L 
converges faster than 1/L & ~ X after determining S. Using 
the calculated values of /(£) and f(L + 2), the exponent 
5 can be estimated by the form 



1 



7(i 



L 



L e-s+x. 



(8) 



The convergence of the size correction is guaranteed by 
the condition Thus the extrapolation of the in- 

dependent exponent 5(L. L + 2) defined by the left hand 
side of (||) gives an estimation of S. 

Note that the method can be easily generalized for the 
behavior around a finite magnetization too, which is de- 
scribed as to — too ~ (H-Hc) 1 / 5 . In this case we have 
only to change the form (J^j into 

f{L) = E(L, M Q + 2) + E(L, M ) - 2E{L, M + 1), (9) 

where M = Lttiq. In addition the method can be applied 
even to gapless cases where H c might be zero. In the 
following argument we don't mention the value of H c 
but we concentrate on the estimation of S. 

For the behavior of the magnetization curve around 
to = 1/2 of the S = 1 antiferromagnetic chain, the in- 
dependent exponent 8{L, L + 2) derived from the form 
(0) using the finite cluster results of E(L, M) up to 
L = 18 calculated by Lanczos algorithm is plotted ver- 
sus 1/(L + 1) in Fig. ^. Fitting the quadratic function 
6{L,L+2) - S + a/(L + l)+b/(L + l) 2 to the data, the ex- 
trapolated value is determined as 6 = 0.99 ± 0.01, based 
on the standard least-square method. The result leads 
to the conclusion 6 — 1 , which is reasonably expected for 
the gapless point m = 1/2. To check the condition (Q) 
E(L, M) — e(m) for m = 1/2 of the 5 = 1 antiferromag- 
netic chain is plotted versus l/L 2 in Fig. ||, where the 
value of c(to) was estimated by fitting of the quadratic 
function of l/L. The plot suggests 9 = 2 in the form (||) 
which is consistent with the prediction of the conformal 
field theory. || Then the condition (0) is satisfied. 
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FIG. 1. L-dependent exponent S(L, L + 2) plotted versus 
1/(Z/+1) for m — 1/2 (solid circle) and m = (solid diamond) 
of the S — 1 antiferromagnetic chain. The extrapolated re- 
sults are 8 = 0.99 ± 0.01 and S = 1.9 ± 0.1, respectively. 
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for (3=2.0 and —0.2 in Fig. ||. The same extrapolation as 
the S = 1 chain results in 6 = 2.03±0.03 and 1.9±0.1 for 
(3 =2.0 and —0.2, respectively. The results are also con- 
sistent with 6 = 2 predicted by some theories discussed 
above. We also have to check the condition ([?]) which is 
8 > 1 because of 6 = 2. In Fig. | E(L, M) - e(m) for 
m = M = Q of the system @ with (3 = 2.0 and -0.2 is 
plotted versus 1/L. It obviously shows a faster conver- 
gence of the size correction for the ground state energy 
per site than which implies that the condition is 

satisfied. 



FIG. 2. Size correction of the ground state energy per 
site E(L, M) — e(m) of the S — 1 antiferromagnetic chain 
plotted versus 1/L 2 for m = 1/2 (solid circle) and m = 
(solid diamond). It converges just like 1/L 2 for m = 1/2, 
while obviously faster than 1/L 2 . 

In Fig. iwe also show the plot of S(L, L + 2) based on 
the form (J|) versus 1/(L + 1) for to = of the S = 1 an- 
tiferromagnetic chain which is more interesting because 
the system has the Haldane gap which vanishes at H c . 
The extrapolated result is 6 = 1.9 ± 0.1 which suggests 
6 = 2, as predicted by the above effective Hamiltonian 
theories. The plot of E(L,M) — e(m) versus 1/L 2 for 
to = in Fig. y obviously shows that the size correction 
in the form (pft decays faster than 1/L 2 in contrast to the 
plot for m = 1/2. It implies 8 > 2 and the condition (0) 
which is 8 > 1 in this case is also satisfied for to = 0. 

Next we investigate the S — 1/2 bond alternating 
chain as another example with the spin gap between the 
singlet ground state and the triplet first excited state. 
The Hamiltonian is defined as 

H = Ho + Hz, 

L L 

H =Y, ^-1 ' S V - /?E ^ • S2 i+*' ( 10 ) 

3 j 
2L 

H z = -H^S*, 

3 

where 2L 5 = 1/2 spins are included in the systems. The 
system has the gap except for /3 = — 1 where it is the 
uniform S =1/2 antiferromagnetic chain. We chose two 
typical values of j3\ (i)/3 = 2.0 and (ii)/3 = —0.2, which 
correspond to the ferromagnetic-antiferromagnetic and 
antiferromagnetic-antiferromagnetic alternating chains, 
respectively. In the latter case particularly the finite size 
effect is larger in the vicinity of the gapless point (3 = — 1. 
Thus we study only for a smaller (3 (=—0.2) than realistic 
cases. The universality argument, however, justifies that 
the critical exponents are independent of (3 except for 
(3 = — 1 , because the system with m = is in a common 
phase for (3 ^ — 1. In order to estimate the exponent <5 
around ttt, = of the system (p^|), the L-dependent expo- 
nent S(L, L + 2) up to L = 1 2 is plotted versus 1/(L+ 1) 
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FIG. 3. L-dependent exponent 5(L, L + 2) plotted ver- 
sus 1/(1/ + 1) for the 5 = 1/2 bond alternating chain with 
13 = 2.0 (solid circle) and (3 = -0.2 (solid diamond). The 
extrapolated results are <5 = 2.03 ± 0.03 and S = 1.9 ± 0.1, 
respectively, //-dependent exponents S± (L, L + 2) associated 
with the magnetization plateau at m = 1/2 are also plotted 
versus 1/(L + 1) for the 5 = 3/2 antiferromagnetic chain with 
the single ion anisotropy D = 8.0; 5+: open circle and <5_: 
open diamond. The extrapolated results are 5+ = 1.98 ±0.04 
and 5_ = 1.99 ± 0.04, respectively. 
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FIG. 4. Size correction of the ground state energy per site 
E(L, M) — e(m) of the S = f antiferromagnetic chain plotted 
versus 1/L for the 5 = 1/2 bond alternating chain (m = 0) 
with (3 = 2.0 (open circle) and (3 = —0.2 (open square), and 
the 5 = 3/2 antiferromagnetic chain (m = 1/2) with the 
single ion anisotropy D = 8.0 (open diamond). It converges 
faster than 1/L in all the cases. We plot the original values 
times 10 4 only for the bond alternating chain with (3 = —0.2. 
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Finally we apply the method to the 5 = 3/2 anti- 
ferromagnetic chain with the single-ion anisotropy. The 
system is described by the Hamiltonian (0) with the 
anisotropy term DJ2j(Sj) 2 added to Ho- Recently an 
argument based on the analogy to the quantum Hall ef- 
fect suggested that the ground state magnetization curve 
possibly had a plateau just at m = 1/2 which corresponds 
to 1/3 of the saturation moment and the singular part of 
the magnetization near the plateau was proportional to 
y/\H — H c \ where H c is the critical field at either edge 
of the plateau. JlTj The plateau was verified to exist for 
D > D c = 0.93 by finite cluster analyses and size scal- 
ing techniques. |18| However, the form of the singularity 
near the edge of the plateau has not been derived by 
any numerical studies on the original Hamiltonian. Thus 
this problem is one of interesting examples to investi- 
gate by the method presented in this paper. We consider 
a sufficiently large D so that the magnetization curve 
has a clear plateau at m = 1/2. The two critical fields 
H c ± are denoted such that the curve has a plateau for 
H c _ < H < H c+ . They can be given by 



c± 



(11) 



although we don't consider the value of H c ± here. To 
investigate the singularity of the magnetization curve, 
the critical exponents S± are defined as 



(12) 
(13) 



To estimate 6± we have only to change f(L) into f±(L) 
defined as f±(L) = ±[E(L, § ±2) +E{L, § ) - 2E(L, | ± 

1) ] and extrapolate the L-dependent exponents S+(L, L+ 

2) defined by the left hand side of the equation (g) using 
f±{L) instead of /(£). In Fig. || we show the plot of 
S±(L,L + 2) versus 1/(L + 1) up to L = 14 for D = 8.0. 
The extrapolated results are 6+ = 1.98 ± 0.04 and <5_ = 
1.99 ± 0.04, which imply S + = <5_ = 2 as suggested by 
the analogy to the quantum Hall effect. We also check 
the condition (Q) by the plot of E(L, M) - e(m) versus 
1/L for M = L/2 and m = 1/2 in Fig. which suggests 
that the size correction decays faster than 1/L. To avoid 
large finite size effects we considered only a large value 
of D (= 8.0) which is not realistic. But it is expected the 
result S — 2 is always true for D > D c (= 0.93) because 
the transition at the critical field belongs to a common 
universality class. 

Recently the magnetization plateau was also investi- 
gated on the S — 1/2 bond alternating chain with the 
next-nearest neighbor interaction Jl{| by the bosoniza- 
tion technique which lead to S = 2 at the edge of the 
plateau at m = 1/4. The result suggests the transi- 
tion belongs to the same universality class as that of the 
anisotropic 5 = 3/2 chain. 



In summary a finite size scaling method to estimate 
the critical exponent 5 associated with the magnetiza- 
tion curve around the critical magnetic field correspond- 
ing to the amplitude of the spin gap of quantum spin 
chains was proposed and applied to the 5 = 1 antiferro- 
magnetic chain and 5 = 1/2 bond alternating chain. In 
addition the behavior of the magnetization curve around 
the edges of the plateau of the anisotropic 5 = 3/2 anti- 
ferromagnetic chain was investigated by the method. All 
the results indicated the same conclusion 6 = 2. 
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